Optimal Control for Electron Shuttling 
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In this paper we apply optimal control technique to derive the control fields that transfer electron 
from one end of a chain of donors or quantum dots to the other. We formulate it as an optimal 
steering problem, and then derive the dynamics of the optimal control. A numerical algorithm is 
developed to efficiently generate control pulses. We apply this technique to transfer an electron 
between sites of a triple quantum dot and an ionized chain of phosphorus dopants in silicon, in 
which we show that a 500 G magnetic field will transfer the electron and two of the four hyperfine 
states of the electron-nuclear spin pair. 



I. PHYSICAL BACKGROUNDS 

The benefits of implementing a quantum computer in 
silicon namely the ability to exploit the techniques 
of the semiconductor industry and long spin coherence 
times, has been offset with challenges including the cou- 
pling of qubits. One mechanism for exchanging quan- 
tum information between qubits is electron shuttling, 
in which spin or charge qubits are literally moved be- 
tween local sites . For dopant spin qubits in silicon Q , 
electron shuttling has been proposed using voltage gates 
and pulses designed analogously to the STIRAP proce- 
dure Q ; this procedure is referred to as coherent tunnel- 
ing by adiabatic passage (CTAP) 0-0|- Similar mecha- 
nisms have been suggested Q for quantum dots in silicon, 
which have also been proposed as qubits (9Ml4|. CTAP 
and other adiabatic procedures avoid populating unde- 
sired sites at any point during the transfer, thereby elim- 
inating issues of decoherence associated with a specific 
site. However, if the source of decoherence is not site- 
specific or is controllable by alternative means, it may 
be useful to approach the state transfer problem for sili- 
con qubits using optimal control theory instead. Optimal 
controls which minimize transfer time or (as will be ex- 
plored here) minimize pulse fluence have been shown to 
minimize decoherence due to additive and multiplicative 
white noise, respectively [H|,[l6j]. Additionally, for qubits 
such as charge qubits in which site-specific decoherence is 
not the main problem, we can use optimal controls in or- 
der to minimize the energy required for gate operations. 

In this paper we will investigate the shuttling of elec- 
trons between the ends of a qubit chain using optimal 
control theory. Depending on whether the chain repre- 
sents dopant spin qubits [f| 0] or lateral quantum dots [1] , 
the control fields will affect the tunnel couplings between 
dopants or the on-site energy of a quantum dot, respec- 
tively. The task is to design some appropriate control 
fields to transport the electron to the end of the qubit 
chain. During this process, quantum information can be 
passed through along the array so as to realize the de- 
sired quantum information processing in the solid-state 



quantum bits. 

This problem can be formulated as an optimal steering 
problem in control theory. The state space consists of all 
the density matrices, i.e., Hermitian matrices with unit 
trace. The system dynamics is governed by Liouville-von 
Neumann equation. The Hamiltonian usually contains 
a drift term, which is a fixed term determined by the 
physical nature of the system, and several control terms 
that can be altered externally as a function of time. The 
objective of steering is to find control fields to transfer 
the system from an initial state to a final state at a fi- 
nite terminal time. This is also known as constructive 
controllability problem [TtJ ■ 

One way to solve this problem is to impose a cost func- 
tion, which may encode minimum energy, shortest time, 
or tracking error penalty. We can then employ standard 
optimal control technique such as Pontryagin Maximum 
Principle (l8j to derive optimality conditions for opti- 
mal control fields. These conditions result in a two-point 
boundary-value problem, which is often difficult to solve 
numerically. 

In the current paper we derive the underlying dynam- 
ics that governs the time evolution of optimal control. 
We take an intuitive approach, which yields the same 
results as from sophisticated Lie-Poisson reduction theo- 
rem [X9j | . For a given initial condition for this dynamics, 
it completely determines the time evolution trajectories 
of controls. Hence, to solve for control fields that achieve 
the desired state transfer, we just need to find an appro- 
priate initial condition. It thus becomes an optimization 
problem on finite dimensional real space. 

To solve the Liouville-von Neumann equation numeri- 
cally, we divide the total time duration into a number of 
small steps and then use piecewise constant functions to 
approximate time-varying control fields. The fidelity of 
the achieved states thus depends on all of the piecewise 
constant control values, which consequently are depen- 
dent on the initial conditions of the dynamics discussed 
above. Using the chain rule, we can obtain the gradi- 
ent of the fidelity with respect to the initial condition 
in closed-form solution. With this approach we can now 
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implement mature gradient type of algorithms to solve 
for control fields that accomplish the steering task. 

To exemplify this approach, we investigate the elec- 
tron shuttling problem for three-donor systems in this 
paper. Our control algorithm derivation can be read- 
ily extended to both systems with more donors and to 
alternative objective functions such as minimizing spin 
decoherence during spatial transfer. We demonstrate the 
efficacy of our control algorithm by applying it to two 
physical systems in Refs. ||| and [|| and then generating 
the control functions to realize the electron shuttling. 



II. MATHEMATICAL BACKGROUND AND 
FORMULATION 

In this section we introduce some necessary mathemat- 
ical background and a general description of the electron 
shuttling problem. 

We consider physical devices with three qubits. This 
type of system has been previously studied in Refs. [f| 
@, [f| . The system Hamiltonian is defined on Lie algebra 
su(3), i.e., all the 3x3 skew-Hermitian matrices. Define 
a basis for su(3) as 



We also use the usual matrix inner product 
(X,Y) = Tr(XY r ). 



(5) 



The dynamics of donor system is determined by the 
Liouvillc-von Neumann equation: 
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The Hamiltonian H of three donors system can be writ- 
ten in a general form as 
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iH = iH + J2 UlXl = aiXi + ^mXu (7) 
i=p+i 



1=1 



1=1 



where Ho is the drift term, and ui are control fields. The 
implementation of the electron shuttling from one end of 
the donor chain to the other amounts to design control 
functions ui to transfer the density matrix p from the 
initial state 
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Define the numbers Cfi- € C such that 
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where K is the dimension of the Lie algebra and K = 8 
for 5u(3). These numbers are the structure constants of 
the Lie algebra su(3) with respect to the basis {Xk}^ =1 . 
Note that the structure constants are antisymmetric in 
all the indices, i.e., 
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for all i, j, k £ {1, . . . , K}. Up to antisymmetry, the 
nonzero structure constants are 
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at time t = T . 



III. OPTIMAL CONTROL AND NUMERICAL 
ALGORITHM 

To solve the state transfer problem we impose 

a cost function to formulate it as an optimal control prob- 
lem, and then use optimal control theory to derive the 
optimality conditions. Based on these conditions, we de- 
velop a numerical algorithm to solve control fields effi- 
ciently. 



A. Optimal control formulation 

For the steering problem ©-I®, in addition to realize 
the desired state transfer, we seek the control fields that 
minimize the following energy type of cost function: 



min J — / L(u)dt. 
«(■) Jo 



(10) 



We thus obtain a standard optimal control problem. To 
solve it, we construct the control Hamiltonian as fol- 
lows Hi: 



H = (9,[-iH,p])+L(v), 



111) 



3 



where \I/ is the costate matrix whose dynamics is given 
by 



* = -^ = -[*#,*]. 



Let 



<f n = {%[X hP ]). 
The control Hamiltonian can then be written as 

K p 



(12) 
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H = -(*, 



E aiXi + y^uiXi, 



l=p+l 



1=1 
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E qf<^f - y~]v,i<i>i +l(u). 



i= P +i i=i 

The optimality condition 



yields that 
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for I = 1, . . . , p. Solving the algebraic equation Eq. (|16|). 
we can obtain the optimal controls as functions of </>/. 
For the sake of completeness, we list all the optimality 
conditions as follows: 

p = ~[iH,p], 

<j = -[iH, *], 

"l 0] [0 0] 

p = 000 , p T = 000 , ( 17 ) 

o oj [o 1 

cj) - dL 
dui 

This is a two-point boundary-value problem, which is 
often difficult to solve numerically. We consider to derive 
the dynamics for 4>i and then transform it into an opti- 
mization problem on finite dimensional space. Taking 
the derivative of <f>i: 

<j>t = (%[X t ,p]) + (%[X t ,p\) 

= ([-iH, *], [X,,p]) + [X u [-iH,p]]) 

= (*,[[*2r,jr,],p]>, (is) 

where the last equality follows from the Jacobi identity. 
For the Hamiltonian given by Eq. ([7]). we have 
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Since the controls can be determined from Eq. (|16p 
as functions of </>, we can further rewrite Eq. (|19[) as a 
compact vector differential equation: 



= ^W, 



(20) 



where = [<fii . . . 4>k\ ■ Note that the dynamics (TH)]) 
is indeed a direct result from Lie-Poisson reduction the- 
orem. See Chap. 13 of Ref. [if 
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and Ref. 

Recalling that the structure constants C* ; are antisym- 
metric in all the indices, we obtain that 
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Hence, 



const, 



(21) 



that is, || 0|| is a conserved quantity along the optimal 
trajectory. 

With the dynamics of cj> in Eq. (|19p . we do not have 
to solve the difficult two-point boundary-value problem 
Eq. (|17p . Instead, we can assume an initial state of 4> 
and solve the differential equation Eq. (TH)]) , If the re- 
sulting control fields ui determined by solving Eq. (|16[) 
arc optimal, they can transfer the density matrix from 
the initial condition to the desired terminal state; other- 
wise, we can calculate the fidelity of the actually achieved 
terminal state p(T) as 



J = Trp T p(T). 



(22) 



Now the state transfer problem can be cast into an 
optimization problem that maximizes the fidelity J in 
Eq. (|22[) by varying the initial condition of <p. 



B. Numerical algorithm 

We can apply a gradient algorithm to find the opti- 
mal initial condition (f>*(0) that maximizes the fidelity. 
To this end, we need to derive dJ/d<j>i(0), the gradi- 
ent of the fidelity J with respect to the initial condition 
4>i(0). Consider piecewise constant control functions. For 
a given time duration [0, T], divide it into N equal inter- 
vals {[tfe^fc+illf^o 1 of length At = t k+1 -t k = T/N, 
where to = and tjy = T. Within the interval [tfe, tfe+i], 
the control field <f>i assumes a constant value 4>i(k). 

By the chain rule, we have 
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We thus need to derive the three differentials in the right 

du m (k) 
d<j>s(k) 



hand side of Eq. ([25]) . First consider ^Prrv- From 



Eq. (jT6j) . we obtain 
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d<j) s \du 



d 
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(24) 



Because the running cost L is defined as a function of 
only the control fields u as in Eq. (fTU)) . is also a 

function of u only. Therefore, we can get by solv- 

ing the algebraic equation Eq. (|24|). For example, when 



L is taken as a quadratic function L = ^(uf + w|), it is 
immediate to obtain = S ms . 

Next consider d f° ffi . Recall that we have written the 
dynamics of ^ in a compact vector differential equation 
(f2"TJ|) . Differentiating both sides of Eq. (|20[) with respect 
to (f>(0), we obtain 



d -F {4 >) = DF { ®- d * 



d<p(0y d(j}(0) KYJ ^'d^O)' 
where the Jacobian matrix DF((f) is given by 
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(26) 



(27) 



From Proposition 6.1 of Chapter 1 in Ref. [22j, we have 
d ■ d d<f> 



d4>{{)) Y dt d4>{{)) ' 



(28) 



that is, it is legitimate to change the order of the differ- 
entials with respect to t and 0(0). Combining Eq. (|25|) 
and Eq. (|2"8|) , we get the differential equation satisfied by 



d 



DF(<j>)- 



dtd4>(0) vvv #(0) 
where the initial condition is given by 

= I. 



d<£(0) 



(29) 



(30) 



Solving the differential equation Eq. (|29|) . we can get 

d<p s (k) 

Lastly derive du dJ ^ ■ From Eq. ([7]), we have 
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iH(k)= aiXt + J2u l (k)X l 
i=p+i 



(31) 



i=i 



and U k = e ~ lH ^ At , where k = 0, . . . , N - 1. Define 
A k = Ul---Ul f _ 1 p T U N - 1 ---U k . 



Then p^r = p(T), Ajy = Pt, and 

J = Tr pxp(T) = TrA N p N = TrA N ^ip N ^i 
= ■ ■ ■ = Tr Aipi = Tr A po- 



(32) 



We need to derive the performance gradient with respect 
to control fields u m {k): 

dJ dTrA k+1 p k+1 d Ti A k+1 U k pkUl 



du m (k) du m (k) 



TrA 



dU k 



k+l 



du m (k) 

From the following formula [1 
d 



p k Ul + U k p k 



du m (k) 
dUl \ 



du m (k) 
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we have 
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dU k 



Al 



e -iH{k)r Xme iH{k)r dT ^ (35) 



du m {k) Jo 
Substituting Eq. (|35|) into (|33|) . we obtain 



dJ 



du m (k) 



TrA fc+1 



Al 



n 



Jo 



= Tr[A fc+1 ,p fe+1 ] / e-^( fc ) T X rne ^( fc )^r. 
Jo 

Since H(k) is a Hermitian matrix, we can diagnolize it 

as 

H(k) = T(k)T(k)T\k), (36) 
where T{k) is a unitary matrix and 

T(k) =diag{ 7l ,7 2j 73}. (37) 

Therefore, 

rAt 

M^X^^dT 



Al 



T{k)e~ tr{k)T THk)X m T{ky r{k) T\k)di 



= T(k) (Tt(k)X m T(k))Q*dTTt(k), (38) 
Jo 
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where denotes the Hadamard product, i.e., element- 
wise product, of two matrices, and = e^ 7b_7a ^ T . For 
7a 7^ 76 ! w e define 



A/ 



s i(7b-7a)At _ ^ 



$ab = / ^abdT 

JO Alb - 7a) 

and for j a = 7b, $ a fc = At. Therefore, 



n 



and 



f/J 



du rn (k) 



T(k)((T\k)X m T(k))Q<P)T^k), 



= Tr([A k+1 , p k+1 ]T(k) 



■((T^(k)X m T(k))(D<£)T\k). (39) 



We have now derived a closed form solution of 
dJ/d(f>i(0), the gradient of the fidelity J with respect to 
the initial condition 0;(O). It is ready to implement gra- 
dient types of algorithms such as the gradient descent al- 
gorithm, the conjugate gradient method, and the BFGS 
method [H]. 

Note that Khaneja et al [25| developed the GRAPE 
algorithm to solve a similar problem. The difference be- 
tween GRAPE and our algorithm is that GRAPE does 
not assume any correlation between neighboring control 
steps and the obtained control fields could thus be oscil- 
latory. As a contrast, our algorithm optimizes over the 
initial condition of a differential equation and therefore 
the resulting control functions have correlation between 
neighboring steps. 

We will then apply our algorithm to two physical sys- 
tems, the triple quantum dot system discussed in Ref. @ 
and the ionized donor chain discussed in Ref. @ . The op- 
timized control fields realizing the desired electron shut- 
tling demonstrate the efficacy of our algorithm. For the 
ionized donor chain we further show that these optimized 
control fields can also provide a high degree of spin state 
transfer when the shuttled electron is coupled to the 
donor nuclei by the hypcrfine interaction. 



IV. TRIPLE QUANTUM DOT 

We now investigate electron shuttling for the triple 
quantum dot system discussed in Ref. [8[- In this sys- 
tem, an electron beginning in the left dot of a three lat- 
eral quantum dot system is moved to the right dot. The 
relative energies of the left and right dots are controlled 
by external gate voltages. The Hamiltonian is given by 



H 



H L {t) Ji 

Ji J-2 

J 2 jUii(t) 



(40) 



where the control fields are the on-site energies and 
fin, and J i and J 2 are the fixed coupling constants 



between nearest neighboring dots. Using the basis in 
Eq. |T]), we can rewrite the Hamiltonian (|4U|) as 

iH = Ji Ai + J 2 A 2 + ~ x 7 H — x & H g ^3 

(41) 

Same as before, we consider the minimum energy cost 
function 

1 f T 

min J = - J (m|(t) + fJ, R {r)) dr. 
The parameters in the Hamiltonian are 

a 1 = J 1 , a 2 = J2 1 Uf=—, u 8 = — -^-j= — . 
Hence 

fx L = 2u 7 , n R =m - V3u 8 , 
and the running cost is 



U8 



The optimality condition Eq. (|16j) becomes 

<p- = 5u 7 — VSug, 08 = 3ug — s/Sur, 
which yields that 



u 7 

and then 



V307 + 8 V&fa + 5</>8 

u 8 = ■ 



fJ-L 



4V3 ' 



12 



2^3 ' m - V3~ 
We can now obtain the dynamics of as 

01 = ^203 - 0407/2 - V3/60408 

02 = -Jl03 - 0508/V3 

03 = -«/201 + Jl4>2 + 0607/2 + V / 3/20608 

04 = 0107/2 + \/3/60i08 - J206 - 2Ji0 7 

05 = 0208/ \/3 + Ji0 6 + J 2 07 - V3J 2 8 

06 = -0307/2 - V3/20308 + J 2 04 - Jl<f>5 

07 = 2Ji0 4 - J 2 5 

08 = \/3J 2 05 



(42) 



The gradient of the fidelity J with respect to 0;(O) can 
be written as 

dJ y-y 1 dJ du m (k) d(f> s (k) 



#;(0) ^ ^ fri du m (k) d0 s (fc) d0 ; (O) 

Y^Vl dJ d <M fc ) 1 dJ d0 8 (fc) 
^ \idu 7 {k) dMO) + ^V3du 7 {k) dMO) 

1 dJ d0 7 (fc) 5 dJ d4> 8 (k)\ 



4V3du 8 (k) d0 ; (O) 12du s (k) d0 ; (O) 
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(A) 



0.4 0.6 
Time(ns) 



(B) 



FIG. 1: (A) State transfer. Blue solid: pn; Green dashed: 
pii\ Red dotted: P33. Ji and J2 were set to -0.1 meV. (B) 
Control functions. Blue solid: [i^; Green dashed: /j,r. The 
pulses are given for iV = 100 segments. (See electronic version 
for color plots). 



The Jacobian matrix DF{4>) is 
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In Fig. [TJ the state transfer in N = 100 segments is plot- 
ted for Ji and J2 set to -0.1 meV and a transfer time of 
1 ns. 
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V. IONIZED DONOR CHAIN 

In this section we apply our control algorithm to the 
ionized donor chain studied in Ref. [f|. The system con- 
sists of three singly ionized phosphorus donors in silicon, 
and one electron shared in the system. The electron be- 
gins on the first phosphorus, site 1, and the pulses are 
designed to move this electron to site 3. The Hamilto- 
nian is given by (setting h = 1): 



H 



-n u (t) 
-n 12 (t) A -Q 23 (t) 
-n 23 {t) 



(43) 



Here the control terms are Q\ 2 and ^23, which are the 
coherent tunneling rate between adjacent dopants. Un- 
der the basis in Eq. ([1} , this Hamiltonian can be written 
as 



iH = -n i2 Xt - tt 23 X 2 



A 



-x 7 



A 

2^3 



A s 



A 



h. (44) 



We can drop the term —I3 as it commutes with all the 
other terms and thus contributes only a global phase. 
Consider the following minimum energy cost function 



min J 



(n 2 12 ( T ) + nl 3 (T))dT 



with the initial and terminal states given in Eqs. p])- 
(|9]). Following the procedure in Sec. IIII1 we have the 
parameters in the Hamiltonian are 



it! = -O12, u 2 = -0 23 , 
and the running cost is 



<i 7 



A 

"2 ' 



as 



2^3' 



T / \ "1 "2 

Liu) = — + — . 
K ' 2 2 

The optimality condition Eq. dTJ5J) yields that 

<pi =ui, 4> 2 =u 2 . 
Hence the optimal controls are given by 

Ol2 = -01, ^23 = -</>2- 



(45) 



For the Hamiltonian ([4"4"]l . the dynamics of in Eq. (fH))) 

becomes 



A 



i = l 



= E 



2^21 



Cfc + ^CS,]^ (46) 



Substituting the values of structure constants in 
Eq. Q into Eq. (|46p . we obtain the complete dynam- 
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ics of i 



02 03 + A0 4 
-0103 - A0 5 



-A0l - 0206 - 20107 

A02 + 0106 + 02 07 - V30208 

02 04 - 0105 

20104 - 0205 

^0205- 



(47) 



To implement the numerical algorithm, it is easy to get 
that 



dJ 
#;(0) 



N-l 

E 

fe=0 



dJ d<j>i{k) dJ dcj> 2 {k) 



du x {k) d0 ; (O) du 2 {k) d0 ; (O) 



(48) 



The Jacobian matrix DF((f>) can be derived from Eq. (|4T|) 

as 
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For the particular Hamiltonian (|44p . we can derive an 
analytic solution for the decomposition in Eq. Q36p. The 
eigenvalues of H(k) are 



A 

7l = -p 
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and the unitary matrix T(k) is 



-Cl 23 /g 2 n 12 /y/gi(gi+A)/2 n 12 /^ 9l { gi - A)/ 2 

-v/(ff i+ A)/(2gT) y/(g i - A)/(2 ffl ) 
^12/32 23 /V3i(5i + A)/2 n 2 3/Vffi(ffi - A)/2) 



Q 2 



where = ^/A 2 + ifi^ + 4fi 2 2 an d 92 = v^la 

In Fig. [3 the optimized pulses in TV = 100 segments are 
shown for A = 2.7 meV, for a transfer time of 1 ns. This 
is in line with the values calculated using tight binding 
theory Q. In Fig. 12(A), the populations of each site 
are shown as a function of time, the populations is fully 
transferred from the first to the third site, and the middle 
site maintains a relatively low populations throughout. 

Since one of the main qubits of interest is phosphorus- 
doped silicon, in which quantum information is stored in 
the spin of the electrons and nuclei, we also used these 
pulses to attempt to transfer a hyperfine spin state along 
with spatially transferring the electron. For these calcu- 
lations, a hyperfine interaction with a splitting of 117.5 
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FIG. 2: (A) State transfer. Blue solid: pn; Green dashed: 
P22; Red dotted: P33. The value of A was set to 2.7 meV. 
(B) Control functions. Blue solid: f2i2; Green dashed: S7 2 3 ■ 
The pulses are given for iV = 100 segments. (See electronic 
version for color plots). 



MHz was used [2g|. The results are shown in Fig. [3] for 
no external magnetic field (A) and a field of 500 G (B). 
The spins of the nuclei at sites 2 and 3 are initialized into 
the j" state, while on site 1 the electron-nuclear system 
is initialized into one of four hyperfine cigenstatcs. The 
distance measure (D) shown in Fig- [3] is a measure of the 
transfer of this hyperfine state, 



D = 1 - \\pt - PhflU, 



(49) 



where ph / is the density matrix for one of the hyperfine 
pure states, and pr is the reduced density matrix of the 
site 3 nuclear and electron spin at the end of the spa- 
tial transfer. The norm used in Eq. (|49p is the induced 
2-norm of the difference matrix, also known as the spec- 
tral norm, which is the maximum singular value of the 
matrix [23] . We have also calculated the fidelity [2^, [2t| , 
the trace distance 12811 . and the Frobenius norm of the 
difference matrix [27j . and while all norms show a sim- 
ilar picture regarding which states are transferred, the 
measure D has the pictorial advantage of following the 
population on the third site when full transfer is occur- 
ring as well as remaining zero when the fidelity is zero, 
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FIG. 3: (A) Transfer of hyperfine eigenstates are attempted 
using the pulses in Fig. [2] The population of sites 1 (solid 
blue) and 3 (dashed green) are shown as well as the measure D 
of the hyperfine state transfer (Eq. (|49[) ) at site 3 (dotted red). 
(B) The same as (A), but with the addition of a magnetic 
field (500 G). The hyperfine states which align the nuclear 
spin with the magnetic field can be transferred, while the 
spin flipped states are energetically forbidden. (See electronic 
version for color plots). 



unlike the Frobcnius norm. At all magnetic fields, the 
| frf) state can be transferred completely from site 1 to 
site 3, because with all of the nuclear spins up the elec- 
trons remain in the hyperfine eigenstate no matter which 
spatial site it is on. Conversely, the hyperfine state | JJ-i) 
cannot be transferred, because this would require flip- 
ping the spins of the nuclei on sites 2 and 3, which is not 
allowed energetically. The electrons which start in this 
state can be transferred spatially, however. Calculations 
have shown that for larger hyperfine constants the spatial 
transfer of the electron can be reduced in the presence of 
the hyperfine interaction. For the linear combinations of 
| ■ftvj.) and | JJ/f), partial transfers can be accomplished 



corresponding to the contribution from the part of the 
linear combination which has the nuclear spin up. This 
can be seen because as the magnetic field is turned on 
(Fig. EJB)) and the linear combinations goes from equal 
to being largely one component or the other, the transfer 
of the component which is mostly nuclear spin up can be 
achieved while transferring the mostly nuclear spin down 
cannot. It should be noted that in the case of a 500 
G magnetic field, two hyperfine states can be transferred 
spatially, suggesting their potential use as a mobile qubit. 



VI. CONCLUSION 

In this paper we have formulated a general problem in 
solid-state electron shuttling as a state transfer problem 
in optimal control theory. We derived the underlying 
dynamical equations that govern the time evolution of 
optimal control fields. This leads to an efficient numerical 
algorithm with a small number of optimizing variables. 
We demonstrate the efficiency of our algorithm by two 
physical examples. 

First, we determined the control pulses for state trans- 
fer between left and right quantum dots in a triple quan- 
tum dot system. Since hyperfine interactions in lateral 
quantum dots can be small, such spatial transfer allows 
the ability to transmit quantum information and possi- 
bly use this ability to couple qubits. Second, we looked 
at the system of an ionized phosphorus donor chain in 
silicon. We again determined optimal control pulses for 
spatial transfer, and then we looked at the transfer of a 
spin qubit using these pulses with hyperfine interactions 
considered. In the 500 G magnetic field case, two hyper- 
fine states of the electron and nucleus can be transferred 
in such a donor chain. As the field is decreased, how- 
ever, only one of four hyperfine states can be spatially 
transferred. Therefore, in order to transfer spin quan- 
tum information in a donor chain in a low-field environ- 
ment, control pulses which take the spin interactions into 
account are desired. 
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